EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP-2001-058 
31 July 2001 



Measurement of the Branching Ratio 
for the Process b — ► r — Z7 r X 

The OPAL Collaboration 



Abstract 

The inclusive branching ratio for the process b — > t~z7 t X has been measured 
using hadronic Z decays collected by the OPAL experiment at LEP in the years 
1992-2000. The result is: 

BR(b -> r~V T X) = (2.78 ± 0.18 ± 0.51)% 

This measurement is consistent with the Standard Model expectation and puts 
a constraint of 

tan f3/M± < 0.53 GeV" 1 
at the 95% confidence level on Type II Two Higgs Doublet Models. 
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1 Introduction 



This paper describes a measurement of the inclusive branching ratio BR(b — > 
t~V t X ), using data taken with the OPAL detector at LEP in the years 1992- 
2000 at e + e~ centre-of-mass energies around the Z resonance. Similar measure- 
ments of BR(b — > r~z7 T X ) have been published previously by the other LEP 
experiments [1-3]. The measurements can be directly compared to the Standard 
Model expectation calculated in the framework of Heavy Quark Effective Theory 
(HQET) dJ, and so can be used to constrain basic parameters of HQET ||. 

A measurement of BR(b — > t~77 t X) is also a probe for the presence of a new 
charged boson coupling to mass. This coupling would increase the branching 
ratio BR(b — > r~T7 T X.) || [7|]. Since a charged Higgs boson contributes at tree 
level, its contribution cannot be easily cancelled by other new particles. This 
can be used to set a limit on a contribution of the charged Higgs exchange. In 
the Minimal Supersymmetric Standard Model (MSSM), however, a region of the 
parameter space is found where one-loop SUSY-QCD corrections could weaken 
the bound |j. 

Many extensions of the Standard Model, like the MSSM, include Type II 
Two Higgs Doublets, where one of two Higgs doublets couples only to down-type 
quarks and the other only to up-type quarks. In these models, tan (3 is the ratio 
of the vacuum expectation values of the two Higgs doublets and is the mass 
of the charged Higgs boson. The decay rate for b — > r~T7 T X can be calculated as 
a function of r = tan/5/M^. A term proportional to r 2 is added to the Standard 
Model decay rate of BR(b -> t~V t X) = (2.36±0.17)%. A value of r = 0.5 GeV" 1 , 
for example, yields BR(b — > r~z7 T X) = (3.61±0.36)% @. The errors are fully cor- 
related and include the experimental uncertainty on BR(b — > e~77 e X) = (10.86 ± 
0.35)% [[| and the theoretical uncertainties. 

Since direct reconstruction of the r lepton in a multihadronic event is not 
possible, other properties of the signal decays have to be exploited. Each event 
is divided into two hemispheres. Hemispheres containing b — > t~77 t X decays 
are characterised by large missing energy, due to the presence of at least two 
neutrinos in the final state. The reconstruction of the missing energy uses the 
OPAL calorimeters and tracking detectors and it relies on the hermeticity of 
the detector. To select a sample enriched in b decays, a b tagging algorithm is 
applied. The b tagging is done in the hemisphere opposite to the signal to reduce 
the dependence on the Monte Carlo simulation of the signal. 

Semileptonic b decays like b — > £~T>X (£ = e, /i) are an important background. 
They are suppressed by rejecting hemispheres with an identified electron or muon. 
The same lepton veto suppresses the leptonic decays of the r lepton in signal 
events, thus selecting mostly hadronic r decays. 
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2 Detector, data set and Monte Carlo samples 



The details of the construction and performance of the OPAL detector are de- 
scribed elsewhere [PDfl . Here only the main components relevant for this analysis 
are described. 

Tracking of charged particles is performed by a central detector, consisting of 
a silicon microvertex detector, a vertex chamber, a jet chamber and ^-chambers fl 
The central detector is inside a solenoid, which provides a uniform axial magnetic 
field of 0.435 T. The silicon microvertex detector consists of two layers of silicon 
strip detectors; for most of the data used in this paper, the inner layer covered 
a polar angle range of | cos#| < 0.83 and the outer layer covered | cos#| < 0.77, 
with an extended coverage for the data taken after the year 1996. The vertex 
chamber is a precision drift chamber which covers the range | cos#| < 0.95. The 
jet chamber is a large- volume drift chamber, 4.0 m long and 3.7 m in diameter, 
providing both tracking and ionisation energy loss (dE/dx) information. The 
z-chambers provide a precise measurement of the z-coordinate of tracks as they 
leave the jet chamber in the range | cos#| < 0.72. 

Immediately outside the tracking volume is the solenoid and a time-of-flight 
counter array followed by an electromagnetic shower presampler and a lead-glass 
electromagnetic calorimeter. The return yoke of the magnet lies outside the 
electromagnetic calorimeter and is instrumented with limited streamer chambers. 
It is used as a hadron calorimeter and assists in the reconstruction of muons. The 
outermost part of the detector is made up by layers of muon chambers. 

Hadronic Z decays collected with the OPAL detector at e + e~ centre-of-mass 
energies around the Z resonance are selected using a standard multihadron se- 



lection pT| . To reduce further the small contribution of Z — > r + r decays an 



additional requirement of at least 7 tracks in each event is imposed. With these 



criteria, the selection efficiency for hadronic Z decays is (98.1 ± 0.5)% |L2| with 
a background of (0.11 ± 0.03)%. After hadronic event selection and Z — > t + t~ 
rejection, the resulting data sample collected in the years 1992-2000 after the in- 
stallation of the silicon microvertex detector consists of 3.70xl0 6 events. About 
11% of the data used were recorded in the years 1996-2000, for calibration pur- 
poses. 

A Monte Carlo sample of hadronic Z decays of about seven times the size 
of the recorded data sample for b flavour events and about the same size as 
the recorded data for the other flavours is used in the analysis. The Monte Carlo 
events are generated using JETSET 7.4 |13[ with the b and c quark fragmentation 



modelled according to the parameterisation of Peterson et al. [H. A global fit to 



OPAL data has been performed to optimise the JETSET parameters [15]. The 



decay b — > r z/ r X is modelled in JETSET using matrix elements neglecting the 



*A right handed coordinate system is used, with positive z along the e~ beam direction and 
x pointing towards the centre of the LEP ring. The polar and azimuthal angles are denoted by 
9 and (j>, and the origin is taken to be the centre of the detector. 
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mass of the final state particles. The energy distribution of the r lepton in the b 
rest frame is therefore reweighted to reproduce the energy distribution calculated 



by including mass effects 0, 16]. The polarisation of the r leptons is simulated 



by interfacing JETSET to the TAUOLA [T7j r decay simulation package. The r 
polarisation is calculated according to j4| with the HQET parameters Ai and A2 
set to zero, corresponding to the free quark model. 

All events have been processed using a full simulation of the OPAL detector 



IS] and the same reconstruction algorithms that were applied to the data. 



3 Event and hemisphere selection 

Each event is divided into two hemispheres using the thrust variable, T, which is 
defined by |TJ 

T -MW)> (1) 

where pi are the momentum vectors of the particles in an event. The thrust 
axis ht is the direction n which maximises the expression in parenthesis. A 
plane through the origin and perpendicular to ht divides the event into two 
hemispheres. 

Signal hemispheres are characterised by large missing energy. Particles es- 
caping close to the beam pipe can fake the large missing energy signature of 
neutrinos. By requiring the polar angle 9 of the thrust axis ht of the event to 
satisfy | cos#| < 0.8, only events well contained in the central detector are ac- 
cepted. Events with a two-jet topology are expected to have thrust values close 
to one and are therefore selected by requiring the thrust T > 0.85. 

The b-tagging algorithm gives a likelihood that a hemisphere originates 
from a b decay by combining the information from a secondary vertex neural 
network, a jetshape neural network and a prompt lepton finder. A cut on the b 
likelihood was applied giving a selection efficiency of 47% and a purity of 92% 
for b flavour hemispheres in the selected sample according to Monte Carlo. 

The hemisphere opposite to the signal hemisphere is used for b-tagging. This 
reduces the dependence on the Monte Carlo description for the signal since the b 
tagging efficiency can be measured directly from the data using events with one 
and two tagged hemispheres. 

4 Missing energy distribution 

In each hemisphere the missing energy E^™ is calculated by 

-^miss - -^boam + -&corr ~ -& v is ■ \ Z ) 
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The sum of the beam energy Eb eSLm and the correction term E corr is the predicted 
energy in the hemisphere. The missing energy E^^ 1 is obtained by subtracting 
the visible energy E^™ 1 . The correction term E corr is determined by exploiting 
the overall energy and momentum conservation in the Z decay. Assuming a decay 
of the Z boson into two particles, the correction term is 

M 2 — M 2 

zp hemi ohcmi /q\ 

where Mh em i and M h em i are the measured invariant masses of the signal hemi- 
sphere and of the opposite hemisphere. 

The visible energy E^?™ 1 in the hemisphere is obtained by summing separately 
the energies of charged particles reconstructed in the central detector and neu- 
tral particles depositing energy in the electromagnetic calorimeter. A matching 



algorithm |21| is used to associate tracks in the central detector with clusters in 
the electromagnetic calorimeters. This algorithm corrects the measured electro- 
magnetic calorimeter cluster energy if a track is matched to the cluster. Tracks 
are counted only if they have been reconstructed using at least 20 jet chamber 
hits, have a p t of at least 120 MeV with respect to the beam axis and a distance 
of closest approach do to the beam axis of less than 2.5 cm. Electromagnetic 
clusters are counted if they have at least 100 MeV energy in the barrel region or 
250 MeV in the forward region of the detector. 

Only hemispheres with missing energy E~ig£- > 5 GeV are used. The cut has 
been optimised to maximise the statistical sensitivity of the fit whilst minimising 
the systematic uncertainty due to the Monte Carlo description of the missing 
energy distribution. The optimisation is discussed in Section ||. 

To remove background from non-r semileptonic b decays, hemispheres with 
identified electrons or muons are rejected. Neural networks are used to identify 
electrons [22|| and muons with momenta larger than 2 GeV. In the selected 



Monte Carlo sample 91% of the hemispheres with prompt electrons or muons 
from semileptonic b and c decays are rejected with this requirement. Tracks 
tagged as originating from photon conversions [^3] in the tracking detector are 
not accepted as electron candidates. 

A sample of 90587 hemispheres remains in the data after applying all cuts. 
Using the Monte Carlo it is estimated to contain about 9% signal events, 3% 
D~ — > t~V t decays, 35% non-r semileptonic b and c decays, 4% light quark 
decays and 49% hadronic b and c decays. A large fraction of the semileptonic 
b and c decays are decays with leptons with momenta below 2 GeV. Using the 
-^mi™ distributions predicted by the Monte Carlo simulation for signal and the 
sum of all backgrounds, a binned maximum likelihood fit is used to determine 
BR(b — > t~z7 t X ) with the fraction of signal events as the only free parameter. 
The result of the fit is 

BR(b -> t-V t X) = (2.78 ± 0.18)%, (4) 
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Figure 1: Distribution of the missing energy in a hemisphere, E^?^, f° r au se ~ 
lected data events and for the Monte Carlo simulation shown on a logarithmic 
scale and on a linear scale. The backgrounds from semileptonic heavy quark de- 
cays into electrons and muons, and from D~ — ► t~V t decays are shown separately. 
The number of b — > r~z7 T X signal events is set to the result of the fit. 
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where the uncertainty takes into account the limited number of Monte Carlo 
and data events [p5fl . The result is obtained by fitting the data samples of the 
different years with the corresponding Monte Carlo samples. The minimum of 
the fits correspond to y 2 values for 34 degrees of freedom in the range 31 to 45. 
The E 1 ^ ™' distribution for the total data sample is shown in Fig. [1|. 



5 Systematic uncertainties 

The measurement of BR(b — > r~v T X ) relies on the Monte Carlo modelling of 
the missing energy distribution for the signal and background events. Other 
systematic uncertainties arise from the reproduction of the selection and veto 
efficiencies for the data in the Monte Carlo simulation and from limited knowledge 
of the branching ratios for decays involving a heavy quark and a lepton. The 
sources of systematic uncertainties are described below and their contributions 
to the total systematic uncertainty are summarised in Table p]. 



Source 


ABR(b -> t~V t X ) 


leptonic E 1 ^™ 1 description 


-0.16% 


hadronic i?^™ 1 description 


+0.15% 


tracking resolution 


±0.04% 


calorimeter resolution 


±0.06% 


b tagging efficiency 


±0.04% 


e veto 


±0.11% 


/i^ veto 


±0.07% 


(xb) = 0.702 ± 0.008 \W[ 


±0.33% 


(x c ) = 0.484 ± 0.008 26 


±0.04% 


b — > r~I7 T X decay modelling 


±0.03% 


semileptonic b decay models 


±0.26% 


BR(b -> £~ z7X) = (10.73 ± 0.18)% [| 


±0.08% 


BR(b -^corc^M ) = (9.69 ± 0.51)% J26| 


±0.05% 


BR(D S - -> r-77 r ) = (7.2 ± 2.3)% ^ g8| 


±0.13% 


/(b^D-) = (18±5)%) m 


±0.10% 


Total systematic uncertainty 


±0.51% 



Table 1: The contributions to the systematic uncertainty on BR(b — ► r u T X ). 



An imperfect modelling of the signal region of the -Eml™ distribution by the 
Monte Carlo simulation can bias the result. We have therefore studied the mod- 
elling of the -EmS 1 distribution using three different signal-depleted control sam- 
ples: a sample enriched in semileptonic b decays, a sample enriched in hadronic 
b decays and a light quarks control sample. In these samples, the ratio of the 



9 



-^mi™ distributions of data and Monte Carlo has been studied. The missing en- 
ergy distributions for the control samples together with the corresponding Monte 
Carlo predictions are shown in Fig. 0. 

Leptonic control sample: A sample enriched in semileptonic b decays is se- 
lected by using the same requirements for the b-tagging as for the signal 
sample (Section |3|). For the leptonic control sample, hemispheres with elec- 
trons and muons are selected and not rejected. To obtain a pure sample, at 
least one electron or one muon candidate with a high value for the output 
of the lepton identification neural networks is required. According to the 
Monte Carlo about 85% of the control sample are semileptonic b and c 
decays. 

The -Eml™ 1 distribution for this class of events is well described by the Monte 
Carlo simulation (Fig. 0a). To estimate the systematic uncertainty associ- 
ated with residual differences, the ratio of the data and the Monte Carlo 
distribution (Fig. 0b) is fitted with a straight line separately for every year of 
data taking. The largest slope a obtained from these fits is used to reweight 
the Ej^™ 1 distribution for semileptonic b decays in the Monte Carlo back- 
ground. The events are reweighted using a weight 

w{E^) = 1 + a(E^™ - E^), (5) 

where -E 1 ^™ 1 is the mean of the missing energy distribution. Using a second 
order polynomial fit instead yields a similar systematic uncertainty. 

Hadronic control sample: To study the description of the missing energy dis- 
tribution for hadronic b decays, a sample is selected by applying b-tagging, 
the lepton veto and by requiring at least 10 GeV energy deposit in the 
hadron calorimeter. These requirements enrich the sample in events with 
large hadronic energy. The agreement between data and Monte Carlo for 
the ^miss 1 distribution is reasonable in the signal region (Fig. 0c, d). The 
systematic uncertainty is estimated using the method described above and 
given in Table [I]. It increases if the cut on is reduced (Table §. To 

keep the quadratic sum of the statistical and systematic uncertainty small, 
the cut i^nfiss 1 > 5 GeV is chosen. All other relevant systematic uncertainties 
show only a small dependence on the E^™} cut. 

Light quarks control sample: A further test of the Monte Carlo description 
of the -Em^ 1 distribution is performed for a light quark (uds) flavour sam- 
ple. The sample is selected by exploiting the fact that events with a hadron 
carrying reconstructed momentum between 0.5 and 1.07 of the beam mo- 



mentum mostly originate from uds primary quarks J3U]. The light quark 
tag is applied in the opposite hemisphere. The E^™ 1 distribution for this 
sample is shown in Fig. ^p. Reweighting the light quark background events 
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-^miss LUI 


stat 


syst 


(GeV) 









±0.14% 


±1.00% 


3 


±0.16% 


±0.45% 


5 


±0.18% 


±0.15% 


7 


±0.21% 


±0.09% 


9 


±0.27% 


±0.08% 



Table 2: The statistical and the systematic uncertainty on the measurement of 
BR(b — ► r~77 T X ) due to the Monte Carlo description for events with hadronic b 
and c decays for different fitting ranges. 



according to Eq. || results in a negligible contribution to the total systematic 
uncertainty. 

Detector effects: The description of the visible energy distribution depends 
on a correct reproduction of the detector resolution in the Monte Carlo 
samples. The tracking resolution has been modified by ±10% to evaluate 



the systematic uncertainty p2| . The energy scale of the electromagnetic 
calorimeter (ECAL) is varied according to the small shift between the mean 
of the ECAL energy distributions in data and Monte Carlo found for an 
inclusive sample of hadronic Z decays. The uncertainty obtained by vary- 
ing the resolution by ±10% is found to be even smaller. The estimated 
uncertainties are listed in Table |l[ 

The following systematic uncertainties are related to the modelling of the 
cuts: 

b-tagging efficiency: The b-tagging efficiency is checked by comparing the 
fraction of events with one and two tagged hemispheres in data and Monte 
Carlo. Using the hemisphere correlation from the Monte Carlo simulation 
a b-tagging efficiency of 45.7% is obtained using the data compared to 47% 
using the Monte Carlo, i.e. the Monte Carlo is overestimating the b-tagging 
efficiency by approximately 3%. This result is consistent with similar recent 
studies p0| . The b-tagged events are reweighted according to this change in 



efficiency and the systematic shift is found to be negligible, so no correction 
is applied. In addition, the uncertainty on the rejection efficiencies for the 
other flavours is evaluated by changing the efficiencies by 10%. Since the 
b-tagging is applied in the opposite hemisphere and the contamination of 
non-b flavours in the selected sample is below 10%, the contribution to the 
total systematic uncertainty is small. 

Lepton identification: An important background with missing energy due to 
neutrinos originates from semileptonic decays of b and c quarks to muons 
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Figure 2: Comparison of the distributions of the missing energy in a hemisphere, 
iffifaj 1 , between the data taken in 1994 and the Monte Carlo simulation for the 
three control samples: a,b) leptonic control sample; c,d) hadronic control sample; 
e,f) light quarks control sample (see text). The number of hemispheres in the 
Monte Carlo is normalised to the number of hemispheres in the data. The lines 
indicate the weight function w(E^^) (Eq. for the fitted slopes a and for a = 0, 
i.e. perfect agreement between data and Monte Carlo simulation for the control 
samples. 
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or electrons. In |24| , 22[ , the systematic uncertainties on the efficiencies for 
identifying muons or electrons determined by the Monte Carlo are found to 
be 5% and 4%. This uncertainty is taken into account in the efficiency for 
rejecting hemispheres with electrons or muons. 

The systematic uncertainties due to the other selection cuts are negligible. 
The remaining errors are due to fragmentation parameters and decay rates which 
have to be taken from other measurements. 

Heavy quark fragmentation and decay modelling: The distribution of the 
missing energy depends on the b fragmentation modelling. The effect 
of the uncertainty is determined by reweighting the Monte Carlo events 
to reproduce the experimental uncertainty on the mean energy of the b 
hadrons, (xt>) using Peterson fragmentation and two other fragmentation 
models |32|, |33|]. This variation produces a different energy distribution for 
the b hadron decay products, resulting in the largest systematic uncertainty. 

Since only about 5% of the sample is expected to be Z—> cc events, the 
uncertainty on (x c ) yields only a small contribution to the total systematic 
uncertainty. 

HQET calculations to order l/m£ show that the r polarisation changes by 
4% compared to the free quark decay model used in the simulation. A 
change in polarisation leads to a different r energy spectrum. To estimate 
the effect of non-zero HQET parameters Ai and A2, the polarisation of the 
r leptons in the Monte Carlo simulation is varied by 4% resulting in a 
systematic uncertainty of 0.03%. 

Semileptonic decay models: A large part of the background is composed of 
semileptonic b decays where the leptons have a momentum smaller than 
2 GeV. Both the fraction of these decays and the shape of the missing 

is 



energy distribution depend on their modelling. The ACCMM model [34 



used for the measurement. Using the ISGW |35| model changes the result 
by +0.14% and using the ISGW** @ model changes the result by -0.26%. 
The larger of the two variations is used as an estimate of the systematic 
uncertainty. 

Uncertainties on decay rates: The branching fractions of the Z boson into cc 



and bb pairs in hadronic events, R c and R^, are taken from |26fl . The varia- 
tion of R c and i?b within their uncertainties yields a negligible contribution 
to the systematic uncertainty. The only significant background involving 
r leptons originates from D~ — ► r~V T decays. Averaging a recent mea- 
surement |27j with the measurement in p8 |, we obtain BR(D,7 — ► t~V t ) = 



(7.2 ± 2.3)%. The branching fraction of b quarks into D~ is f(h — > D~) = 
(18 ± 5)%) [p9fl . These uncertainties and the uncertainties on the semilep- 
tonic branching fractions BR(b — > £~VX) and BR(b — > c or c — > £z/X ) have 
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been taken into account by reweighting the Monte Carlo events using the 
branching fractions and the uncertainties given in Table |l|. 



6 Results 

Using about 3.70 million hadronic Z decays we have measured the inclusive 
branching ratio: 

BR(b -> r~z7 r X) = (2.78 ± 0.18 ± 0.51)%. (6) 

The size of the total uncertainty is similar to the uncertainties of the measure- 
ments in [1-3]. 

A contribution from charged Higgs decays is expected to enhance the branch- 
ing ratio compared to that in the Standard Model. Since we have found no large 
enhancement of the branching ratio compared to the Standard Model prediction 
BR(b — ► r~77 T X) = (2.36 ± 0.17)%, we set a constraint on such a contribution for 
Type II Two Higgs Doublet Models 0. 

The t polarisation depends on the Higgs contribution and can be calculated 
as a function of r = tan/3/M H ± 0. This is taken into account by an iterative 
procedure for the limit calculation. A limit of r < 0.52 GeV" 1 is obtained 
assuming r polarisation as for the Standard Model decay. The polarisation of 
the r leptons in the Monte Carlo simulation is changed according to this value 
of r and the limit is recalculated. The resulting limit is 



M H ± 

at 95% confidence level. 



< 0.53 GeV" 1 (7) 
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